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Neutrinos and antineutrinos emitted from a core collapse supernova interact among them- 
selves, giving rise to collective flavor conversion effects that are significant near the neu- 
trinosphere. We develop a formalism to analyze these collective effects in the complete 
three-flavor framework. It naturally generalizes the spin-precession analogy to three fla- 
' vors and is capable of analytically describing phenomena like vacuum/MSW oscillations, 

I synchronized oscillations, bipolar oscillations and spectral split. Using the formalism, we 

demonstrate that the flavor conversions may be "factorized" into two-flavor oscillations with 



(N 



hierarchical frequencies. We explicitly show how the three-flavor solution may be constructed 
by combining two-flavor solutions. For a typical supernova density profile, we identify an 
' approximate separation of regions where distinctly different flavor conversion mechanisms 

operate, and demonstrate the interplay between collective and MSW effects. We pictorialize 



our results in terms of the "ea-eg triangle" diagram, which is a tool that can be used to 
visualize three-neutrino flavor conversions in general, and offers insights into the analysis of 



the collective effects in particular. 
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INTRODUCTION 



Neutrinos emitted from a core collapse supernova carry information about the primary fluxes, 



^ , neutrino masses and mixing, and SN dynamics 



l|, Neutrinos, produced in the region of the 
' neutrinosphere, freestream outwards and pass through the core, mantle and envelope of the star. 

The drastically different environments in these regions, consisting of varying densities of ordinary 
matter, neutrinos and antineutrinos, affect flavor conversions among neutrinos. Multiple shock 
fronts and turbulence generated during the explosion may also affect the neutrino flavor conversions. 

The traditional picture of flavor conversions in a SN is based on the assumption that the effect of 
neutrino-neutrino interactions is small. In this picture, neutrinos that are produced approximately 
as mass eigenstates at very high ambient matter density in the core propagate outwards from the 
neutrinosphere. Flavor conversion proceeds most efficiently at the electron density corresponding 



to the MSW resonance 



I\4 



J|. The outcoming incoherent mixture of vacuum mass eigenstates 



is observed at a detector to be a combination of primary fluxes of the three neutrino flavors. 



This scenario of resonant neutrino conversions in a SN has been studied extensively to probe 
neutrino mixings and SN dynamics. The work has focussed on the determination of mass hierarchy 



and signatures of a non-zero ^13 



pass through matter 



independent signatures 



11 



earth matter effects on the neutrino fluxes when they 



jshock wave effects on observable neutrino spectra and their model 



la ]. Recently, possible interference effects for multiple resonances jl7l |. 



18h2C1I|. and time variation of the signal 



the role of turbulence in washing out shock wave effects 
21[ have also been explored. 

However, neutrino and antineutrino densities near the neutrinosphere are extremely high 
^2q30-35 pgj, gj^^s-j^ which make the neutrino-neutrino interactions significant. Such a dense gas 
of neutrinos and antineutrinos is coupled to itself, making its evolution nonlinear. The flavor off- 

23 LA formalism 



24 



261 ]. where a 



diagonal terms can be sizeable, and significant flavor conversion is possible |22l . 
to study flavor evolution of such dense relativistic neutrino gases was developed in 
set of quantum kinetic equations for their evolution were written down. These equations have been 
studied in detail, though mostly in the two-flavor approximation, and the nature of flavor evolution 



has been identified 



27 



collective flavor conversion, i.e. neutrinos of a . 
oscillations 



301 ]. It was eventually understood that a dense gas of neutrinos displays 

I energies oscillate together, through synchronized 



311 ] and/or bipolar oscillations 



33] ■ Another remarkable effect of these interactions 



is a partial or complete swapping of the energy spectra of two neutrino flavors, called step-wise 
spectral swapping or simply spectral splits, as the neutrinos transit from a region where collective 



3614401]. Recent two-flavor 



41 



effects dominate to a region where neutrino density is low 

The nonlinear effects in the context of SNe were considered in 
simulations showed that the collective effects affect neutrino flavor conversions substantially 
42l |. Different collective flavor transformations seem to play a part in different regions of the 
star. Many features of the results of these simulations can be understood from the "single-angle" 
approximation, i.e. ignoring the dependence of the initial launching angle of neutrinos on the 



evolutions of neutrino trajectories 



32 



35 



4314481]. Angular dependence of flavor evolution can 



give rise to additional ang le dependent features observed in two-flavor simulations 
to decoherence effects 



2g 



41 



il, or 



491 ]. For a realistic asymmetry 



aetween neutrino and antineutrino 



fluxes, such angle dependent effects are likely to be small 



47 



481 ]. Recently collective effects have 



also been numerically investigated in the context of the neutronization-burst phase of 0-Ne-Mg 
supernovae jsol . The impact of these nonlinear interactions has also been studied in the context of 
cosmological neutrino flavor equilibration in the early Universe where the synchronized oscillations 



play a significant part 



51 
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Most of the analytical results in this area are in the two-flavor approximation, where the equa- 
tions describing flavor transformations are similar to the equations of motion of a gyroscope. As 
a result, two-flavor results are now fairly well understood analytically, with the exception of pos- 
sibly two issues, viz. decoherence (or lack of it) for asymmetric systems 47[, and the existence 



of the antineutrino spectral split [48]. In this work we focus on the effects of the mixing of the 
third flavor. We present an analytical framework to study three-flavor effects and demonstrate an 
approximate factorization of the full three- flavor problem into simpler two- flavor problems, when 
AmQ <^ |Am^tn^| and ^13 <^ 1. We also develop the "ea-eg triangle" diagram, a tool that can be 
used to visualize the three neutrino flavor conversions in general, and allows one to gauge the extent 
of additional effects of the third flavor. We numerically study collective flavor transformations for 
a typical SN density profile, identify regions where different flavor conversion mechanisms operate, 
and explain the features of the spectra using our formalism. 

The outline of this paper is as follows. In Sec. [Ill we review the equations of motion of a 
dense gas of neutrinos in steady state. We specialize these equations to spherical geometry in 
the single-angle half-isotropic approximation, and write the three-flavor analogue of the gyroscope 
equations, by introducing the eight-dimensional Bloch vectors. We recover the two-flavor limit of 
those equations, and recognize an approximate factorization of the three- flavor problem (reminis- 
cent of the H-L factorization in the standard picture) into smaller two-flavor problems. We show 
that survival probabilities can be written down in a simple form, purely in terms of the solutions 
to the two- flavor problems, as long as the frequencies governing the oscillations are hierarchically 
separated. In Sec. IIIH we illustrate the above factorization for vacuum/MSW oscillations as well 
as collective synchronized oscillations. We also explain the three- flavor features of bipolar oscilla- 
tions and spectral splits qualitatively and pictorially. In Sec. \IV\ we calculate the flavor conversion 
probabilities numerically for a typical SN matter profile, and identify the additional features due 
to the third neutrino. We conclude in Sec. |V] with a summary of our results and comments about 
directions of future investigation. 



II. THREE-FLAVOR FORMALISM 



In this section we derive the steady state equations of motion for an ensemble of dense gas of 
three- flavor neutrinos, as a straightforward generalization of the corresponding equations in the 
two- flavor case. 
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A. Hamiltonian and Equations of Motion 



We denote a neutrino of momentum p at time t at position r by z^(p, r, t), and work in the 
modified flavor basis {ue,iyx,i^y) defined such that (fg i^x I'y)'^ = -R23(^23)('^e i^fi ^t)"^ ■, where -R23 is 
the rotation matrix and ^23 the neutrino mixing angle in the 2-3 plane. ^ In this basis, the density 
matrix for n,y(p, r, t) neutrinos with momenta between p and p + dp at any position between r 
and r + dr may be written as 

1 



'^|i^(P,r,t))(i/(p,r,i)|c«/; 



(1) 



where a, /3 = e,x,y and the summation is over all nu{p, r, t) neutrinos. Note that the density matrix 
is normalized to have unit trace, but the neutrino density itself is n,^{p,r,t), which typically falls 
off as 1 /r^ from the source. The number density of neutrinos with flavor a is obtained through 

riu^ (p, r, t) = n^{p, r, t)p^^^^ (p, r, t) . (2) 

We define the analogous quantities for antineutrinos similarly, but with a reversed order of flavor 
indices to keep the form of equations of motion identical [2^ , and denote them with a "bar" over 
the corresponding variables for neutrinos. 

The effective Hamiltonian in the modified flavor basis for neutrinos z^(p, r, t) of energy E ^ p = 
IpI in vacuum is 

H.oac{p) = UM^Uy2p , (3) 
where the masses and the mixing matrix are parametrized as 

M = Diag(mi,m2,m3) , (4) 



U = R 



23(^23)-R23(6'23)-Rl3lC'13j-f1.i2l,C7i2j , 



)Rl 



(5) 



with Rij being the rotation matrices in the i-j plane. In this work, we take the value of the 
CP-violating phase in neutrino sector to be zero. Now H^acip) may be explicitly written as 



2p 



( sfg C13S13 ^ 



\ C13S13 Cfg j 

„2 „2 



+ 



2p 



C12C13S12 



C12C13S12 



•-12 



2 \ 

-C13S12S13 * 



-C12S12S13 



\ -Cl3Sf2Sl3 -C12S12S13 ■sf2Sl3 / 



(6) 



^ This basis has also been denoted in the literature as {ve, u^i, u^i) 
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where Am?^- = m| — mf and other symbols have their usual meanings. In matter, neutrinos feel 



the Mikheyev-Smirnov-Wolfenstein (MSW) potential 



due to charged leptons 



V{r, t) = V2Gf n^- (r, t) Diag(l, 0, 0) 



(7) 



that adds to the Hamiltonian, where n^- (r) is the net electron number density at r. The effective 
Hamiltonian also includes the effects of neutrino-neutrino interactions, which to the leading order 



in Gp depend only on forward scattering and contribute 



24 



26| 



(q, r, t)p{q, r, t) - n^{ci, r, t)p{q, r, t) 



(8) 



The interaction strength is dependent on the angular separation of the momenta of the interacting 
particles, and is given by Kpq = 1 — cos6'pq, where 6'pq is the angle between p and q. 
The equation of motion for the density matrix is 



|p(p,r,t) 



H{p,r,t),p{p,r,t) 



+ ^p(p,r,t) 



(9) 



In the steady state (no explicit time dependence in the Hamiltonian and initial conditions) and 
ignoring external forces (terms depending on dp/dt), we can drop the time dependence in the 
problem. Writing v = dr/dt we have the equations of motion for /)(p, r) and p(p, r) as joi] 



V • drp{p, r) 



V • drpip, r) 



+ H^acip) + V{r) + H^^ip, r),p{p, r) 



Hvacip) + V'(r) + H^^{p,r),p{p,r) 



(10) 
(11) 



41 



he effect of ordinary matter can be "rotated away" by working in the interaction picture 



44( 1 . We employ an operator 0(r) under which a matrix A transforms to 

A'^'^r) = 0{r)AO-\r) , 

where 

0(r) = exp (^i dr'Vir')^ . 
This choice simplifies the equations of motion by removing the matter term, giving us 
v9ry9^"*(p,r) = -i +HZi{p,r)+HZ\p,r),p'^\p,r) 

V V'^*(p,r) = -i\-HZi{p,r)+Hll\p,r),r\p,r) 



(12) 

(13) 

(14) 
(15) 



We assume that the density of /i* and is negligible, and that and Ur feel approximately identical potentials, 
which have been taken to be zero by convention. An analysis of collective effects including a /i — r potential has 
recently been carried out [goI ]. 
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The transformation by 0(r) leaves diagonal entries of p(p,r), p{p,r), H^ac{p) and H^^{p,r) un- 
changed, but the off-diagonal entries become r-dependent. For example, if V{r) varies adiabatically 
and only in the radial direction, the vacuum Hamiltonian changes according to Eq. ()12p as 

^2 



~2r 



V{r) 



V{r),H,ac{p) 



+ ... . 



(16) 



We know that V{r) is a diagonal matrix, so only the off-diagonal elements of Hq{p) are affected 
by the transformation. The final observables we are going to be interested in, the number fluxes of 
neutrino flavors, involve only diagonal elements of the density matrix [see Eq. ([2])], so the interaction 
basis is well suited for our purposes. 



B. Spherical Symmetry and Single-angle Equations of Motion 

The interaction term H^i, in Eq. ([8]) depends on 0pq, i.e. the angle between the momenta 
of interacting neutrinos. Thus while performing the angular integrals therein, the dependence of 
the neutrino flux on all angular variables must be taken into account. This makes the problem 
quite complicated, and an approximate treatment is needed in order to gain useful insights. Two 
levels of approximation have been considered in literature, viz. multi-angle and single-angle. In the 
multi-angle approximation, azimuthal symmetry about the axis defined by the source and observer 
is usually assumed, but not complete spherical symmetry. This essentially captures the effects 
of correlations between trajectories with different initial launching angles. The effects of such 
correlations can have interesting consequences which have been explored in detail 



2S 



41 



42 



47 



49l |. In the single-angle approximation, it is assumed that the flavor evolution does not significantly 
depend on any of the angular coordinates (i.e. the evolution is spherically symmetric), and thus we 
can integrate over all angular degrees of freedom trivially. One must then choose a representative 
value for cos^p, which we take to be 1/2. 

We assume half-isotropic emission from a source of radius rg, as defined in and write 

n,,(p,r) = nuip,r) = n^{p,ro) rl/r'^ , (17) 
p(P,r) = p{p,r) . (18) 

In the steady state, the fluxes of neutrinos and antineutrinos can be written as 

= J dp 2V 47rrg n^{p,ro) , (19) 
= J dp 27r/ 47rrg MP,ro) , (20) 
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the total flux being $ = + 

A further "unification" in the notation for neutrinos and antineutrinos is possible by noting 
that their equations of motion, i.e. Eqs. (fTOll and (fTTIl . differ only in the sign of H^acip)- This 
suggests a change of variables from p to 

CO = \Aml3\/{2p) . (21) 



Using the same convention as 3J|, we define for neutrinos 

n„{uj,r) = ny{ p{uj), r) , p{uj,r) = p{ p{uj), r) , (22) 

and for antineutrinos 

n„{-u;,r) = p{uj), r) , p{-uj,r) = p{ p{uj), r) . (23) 

The negative values of u thus correspond to antineutrinos. Then we need to solve only for p{uj, r), 
albeit at the cost of extending the domain of to to both positive and negative values. This simplifies 
the H,yi,{p,r) term in Eq. ([8]) to ^ 



Huuir) = p{r) i du f{u) p{u},r) sgn{uj) . (24) 

J —oo 

in terms of the distribution function 

fiu) = ^\^^^^^n^iu,ro), (25) 
normalized as dco f{Lo) = 1, and the "collective potential" 

P{r) = Po g{r) . (26) 
Here is the collective potential at the neutrinosphere: 

and the "geometric dilution factor" g[r) is given by 

di.'^) = / (i(cos6'q) (1 - cos6'qCos6'p; 



cos 6»p=l/2 



^ Note that _ff^y(p,r) depended on p only through the direction of p. This dependence no longer survives in the 
single- angle approximation. 
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The geometric dilution factor equals unity for r = vq, whereas at large rg, it decreases as 
The decrease of neutrino densities from a finite source accounts for a factor of whereas the 
additional dilution factor of approximately arises from the integral in Eq. (I28[) . due to the 
decreasing angle subtended by the source and reduced collinearity, which are encoded in the limits 
and the integrand respectively 43|]. Note that the exact numerical factors depend on the choice of 



cos 9p. 

The total flux <I> remains conserved as long as there is no explicit temporal variation of the 
overall luminosity. We work in the steady state approximation and assume the luminosity to be 
constant in time. Slow variations in it may be taken into account by including an additional time 
dependent factor. Note that /(w) is independent of r, which embodies the fact that the normalized 
neutrino spectrum does not change. Using Eq. ([2]), we can also write the flavor dependent w-spectra 
fu^{yj,r) as 

US^^r) = f{u)p^^^^{uj,r) . (29) 

Note that fuei^^ ^) contains the spectra of both I'e and Ue, and depends on r only through pi,^i,^{uj, r). 
It would be a constant on the trajectory if there were no flavor evolution of p^^^^{u},r). For later 
use, we define the energy integrated neutrino fluxes for each flavor as 

/•oo 

$,,(r) = ^ / dioUX^^r) , (30) 

= ^uSr) + ^uAr)+^uy{r) , (31) 

cl>p,(r) = / duf,Suj,r) , (32) 



$p = $p,(r) + $p^(r)+cl>p^(r) . (33) 

With these approximations, the problem is reduced to an ordinary one dimensional problem 
along the radial direction. We denote the derivative with respect to r using a "dot", and using 
Eqs. (llOp and (jlip . arrive at the single-angle equations of motion 



VrpiyJ-, r) = —i 



+ H^aciu}, h) + Vir) + H^^{r),p{oo, r) 



(34) 



where h = sgn{Am'f^) encodes the dependence on mass hierarchy. Here 



Vr = \ll-^sm^epiro) (35) 



is the radial velocity of the neutrino. Note that for r ^ ro, we have ~ 1. Since the flavor 



conversions due to collective effects start becoming significant only for r > 4ro 



47 



48l |. for our 
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analytic approximations we take Vr = 1, making Eq. ([3 



p(w,r) 



+ H^ac{^, h) + V{r) + H^^{r),p{uj, r) 



(36) 



We have thus used the spherical symmetry of the problem, and the simple energy dependence, 
to rephrase the equations of motion in a somewhat simpler form. This single-angle approximation 
is probably crude, but it has been shown in numerical simulations (for two flavors) that this 



approximation seems to work reasonably well 



48l |. It also seems that the multi-angle effects are 



suppressed when the neutrino and antineutrino spectra are not identical [47[. We assume the above 
results to hold true for three flavors as well, and ignore multi-angle effects in this work. Thus, for 
an analytical understanding of various flavor conversion phenomena associated with this system, 
we confine our discussion to the steady-state single-angle half-isotropic approximation that we have 
outlined above. 



C. Bloch Vector Notation 



In the single-angle approximation, it is useful to re-express the density matrices and the Hamil- 
tonian as Bloch vectors. The idea, analogous to the two-flavor case, is to express the matrices 
in a basis of hermitian matrices, and to study the motion of the vectors constructed out of the 
expansion coefficients (which are called the Bloch vectors). In our problem, we choose the basis 
consisting of the 3x3 identity matrix /, and the 8 Gell-Mann matrices given by 
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1 
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i 





, A3 







-1 
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Ay = 








—i 
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0-2 











(37) 



which satisfy the SU{3) Lie algebra 



[Aa, Afe] = ifabc Ac , 



(38) 
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where o, b, c are integers from 1 to 8. Note that the normahzation for the matrices is chosen such 
that 

T>(A,A,) = 28ab . (39) 

The structure constants fabc are antisymmetric under exchange of any two indices and are specified 
by 

/l23 = 2 ; /l47) /l65) /246) /257) /345) /376 = 1 hrS, /458 = ■ (40) 

Note that basis of traceless matrices can be expressed as a semi-direct sum of 

K = {Ai,A2,A3,A8} and Q = {A4, A5, Ag, A7} , (41) 

i.e. for Ka G K and Qa G Q we have 

[Ka, Qfe] G K , [Qa, Qb] G K and [Q„ Kf,] £ Q . (42) 

In fact this is not the only choice of K and (Q that has this property. In addition to the decompo- 
sition 

K^" = {Ai,A2,A3,A8} and Q^" = {A4, A5, Ag, Ay} , (43) 

as above, we could also choose 

K"^^ = {A3,A4,A5,A8} and Q^^ = {Ai, A2, Ag, A7} or (44) 
K"?^ = {A3,A6,A7,A8} and Q"^ = {Ai, A2, A4, A5} , (45) 

which satisfy the conditions in Eq. (142p . The meaning of the superscripts {ex, ey, xy) on K and Q 
will become clear later. 

Using the basis matrices / and A^, we now express any 3x3 hermitian matrix X as a vector 
X in the 5C/(3) generator space (with unit vectors ej) as 

X = ^Xo/ + ix-A. (46) 

We call X the Bloch vector corresponding to the matrix X. The vector X must lie completely 
within an eight-dimensional compact volume, called the Bloch ball, whose various two-dimensional 
sections are shown in Fig. [TJ We say that the vector X is contained in K^^' (K^^, K^'^) if the matrix 
X can be expressed solely as a linear combination of A^ G K^^ (K^^, K^^). 



11 



Type I 2/^:1 




2/V3 



Type; II 




Type III ^/V^ 




A, e {Ai,A2,A3}^ A, = As 



A,: e {A-i, A5; Ag, A7}, A; = iAs A^ e {A4, A5, Ae, A7}, A^ = As Other combinations 



FIG. 1: The shape of the Bloch-ball for a vector A^e^ (Figure taken from [62| ) 



We reparametrize our equations using Eq. (|36|) , and define the Bloch vectors corresponding to 
the density matrices as 



(47) 



Note that A is an eight-vector of 3 x 3 matrices. The scalar Po(w,r) and the polarization vector 
P(ti;, r) encode the flavor content of neutrinos of energy -p = | Am^3|/(2ti;) at a distance r for a; > 0. 
The negative values of w encode the same information for antineutrinos. Since p{uj, r) has been 
normalised to have unit trace by definition, Po(cL!,r) is equal to one. We will therefore not worry 
about the zeroth component of the polarization vector henceforth. For a pure state, P(a;,r) lies 
on the boundary of the shaded region in Fig. [H and has the magnitude 2/\/3. For a mixed state, 
the magnitude of P(a;,r) is smaller and the vector lies within the shaded region. 

We assume that all neutrinos are produced as flavor eigenstates, i.e. the primary flux con- 
sists of (p, To) and ?T-jy„(p, tq) with energy p. The initial density matrix p{p,rQ) is therefore 
Diag^n,yg(p, ro),n,y^(p, ro),n,^j^(p, ro)^ , and similarly for antineutrinos. The initial polarization 
vector may be written as 



T,/ ^ fue ^o) - fu, (.(^, ro) . , fu, (w) + fu^ (uj) - 2/^ {uj, ro) , 
P(a;,ro) = jj—^ 63 H ^ , eg . 



(48) 



The polarization vector P(a;,r), when projected onto the 63-65 plane, must lie within the 
triangle in Fig. [2l where we show a representative P(c<j, r) projected on the 63-63 plane. The pure 
electron flavor is represented by 



63 + 



The or z^g content with energy p at position r is given by 

/ ^ _ ny^jp^r) _ fu,{^,r) _ 1 ^ P • 6e _ 4 
Pu.uAP.r)- ^^^^ - ^^^^ -3+ 2 



(49) 



(50) 
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3 



e. 



y 



FIG. 2: The projection of a polarization vector P on the es-eg plane 



The projection of P on is thus related to Pu^UeiPii^) = US'^^'^)/ fi'^) ^ above. The same 
quantity can be easily visualized from the figure as (ie/\/3, where is the distance of the tip of P 
from the side of the triangle that is perpendicular to Cg (as shown in the figure). The number of 

and Uy are also similarly calculated. Negative values of oj encode the same information for the 
antineutrinos. This gives a simple pictorial way to represent the flavor content of the ensemble by 
plotting the tip of the projection of P(aj, r) on the 63-63 plane. ^ 

For the mass term in the Hamiltonian, we have 



Note that oj for neutrinos is always positive in this convention, and the negative sign of Amfg for 
inverted hierarchy is absorbed into B. The terms involving e = Am^g/Amfg arise from the mixing 
of the third flavor, and the three-flavor effects enter through them. The sign of e is positive if 

the mass hierarchy is normal (Amf3 > 0) and negative otherwise. This, along with the overall 
sign due to h, guarantees that the contributions from Amf2 always have the same sign. Note that 
B2, B5, By vanish in the absence of CP- violation. 




(51) 



where 




(52) 



* Note that probability conservation in tliis representation corresponds to tlie tlieorem tliat tlie sum of tlie lengtlis 
of perpendiculars dropped from any point inside an equilateral triangle to tlie three sides is a constant. 
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The MSW potential defined in Eq. may be represented as 

V{r) = AwQloI + ^L-a) , (53) 

where A(r) = \plGp nf.-{r). The vector L parameterizes the effect of background electrons, and 
is given by 

L = eg + eg/Vs . (54) 
The liyy(r^ term defined in Eq. ()24p can also be simply written as 

F,,(r) = /.(r) Q Do / + ^ D(r) • , (55) 

where D(r) is defined as 

D(r) = y '^'^ A^) P(w,r) sgn{u:) . (56) 

In the next section, we shall express the evolution equation, i.e. Eq. ()36p in terms of the Bloch 
vectors P(a;, r), B(a;, /i), L and D(r). 

D. Generalized Gyroscope Equations 

We have expressed our problem in terms of the eight-dimensional Bloch vectors, and now we 
shall see that the equations of motion formally resemble the equations of a gyroscope. To make this 
apparent, we define x as a generalized "cross product" josl with fahc as the structure constants, 
instead of the usual e^fec that appears in the two-flavor approximation, e.g. 

8 

B X P EE ^ /afecBaPf, • (57) 
a,fe=l 

This makes it possible to write the equations of motion, i.e. Eq. p6p . compactly as 

P(w,r) = ( a;B-FA(r)L-F//(r)D(r) ) X P(tj,r) = H(tj, r) x P(a;, r) , (58) 



where P(a;,r), B, L, D(r) are defined in Eqs. ()17|) . (j52]) . ([51]) and ([56]) respectively. The couplings 
w, /u(r) and A(r) are defined in Eqs. (j2T]) . ([28]) and ([53]) respectively. Equation (j58]) resembles the 
equation of a spin in an external magnetic field, or equivalently, that of a gyroscope. We must 
remember that this similarity is purely formal, because unlike in the two-flavor case, we cannot 
write an arbitrary Bloch vector as a linear combination of two Bloch vectors and their cross product. 
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We shall show in Sec. HIE] that under certain approximations these generalized gyrosope equations 
can be given a geometrical interpretation. 

The effects of the matter term A(r)L in Eq. ()58p can be rotated away by going to the interaction 
frame as described in Eq. (|12p . where a matrix A becomes = OAO~^. In order to determine 
the Bloch vector corresponding to A*"*, we equate 

^I + j2^ = 0A0-K (59) 

a=l 

Multiplying both sides by Aq and taking trace, we get 

^mt ^ TiiAaOAO-^) , (60) 
where we have used Tr(AaAfc) = 26ab- In particular, the Bloch vector B*"* may be written using 



Eqs. UM and §0^i as 

B^"*(r) = BicosC(r) ei + BisinC(r) ea + Bs es 

+B4 cos C{r) 64 + B4 sin C{r) + Bg eg + B7 67 + Bs eg , (61) 

where (^{r) = V{r')dr' . In dense matter, B^"*(r) oscillates rapidly with the freq uency ~ 



mimicking a suppression in the relevant mixing angles as in the two- flavor case [32]. 
We also define the "signed" and "unsigned" n*^ moments (with n > 0) of P(a;,r) as 

D(")(r) = j duuj'' f{uj) P(w,r) sgn{uj) , (62) 
S(")(r) = J doju'' f{u;)P{uj,r) . (63) 

Note that D^"-* is same as D, and we will therefore refer to S^*^^ as S. The evolution of these 
moments are governed by 

D(")(r) = B X D("+i)(r) + (^A(r)L + /i(r)D(r)^ x D(")(r) , (64) 



S(")(r) = B X S("+i)(r) + (^A(r)L + /i(r)D(r)^ x S^") 



(r) . (65) 



We see that the higher moments turn up in equations of motion the lower moments. If we take the 
dot product of Eq. §^ with D(")(r), and of Eq. with S(")(r), we get 

a,.|D(")(r)|2 = D(")(r) • B x D("+i)(r) , 

5^|S(")(r)|2 = S(")(r) • B x S("+i)(r) . (66) 

The above dependence of the moments on r implies that there is likely to be a redistribution of 
flavor as a function of w. It will be interesting to investigate if these moment equations can be 
used to predict the nature of the redistribution of flavor spectra. 
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E. Heavy-Light factorization of dynamics 

The three-flavor dynamics in the traditional matter-driven scenario can be factorized into the 
so-called "heavy" and "light" MSW resonances that occur at densities corresponding to Amatm 
and Attiq respectively. Appropriate combination of the effective two-flavor dynamics in these two 
sectors approximates the three-flavor result reasonably well. We now proceed to illustrate a similar 
simplification for collective effects as well Let us first introduce the notion of "heavy" and "light" 
subspaces of the Bloch-sphere. In the K-Q decomposition shown in Eq. ()44p . the vectors contained 
in are termed "heavy" (H) whereas those contained in (Q^^ are termed "light" (L). A general 
vector X may be decomposed as 

X = . (67) 

In particular, B in Eq. ()52p may be expressed as B = B^ -1- B^, with 

hB" = (^sls - e(c\^ - c\^s\^ 63 + (1 - es?2) sin 2^13 64 

+ (^(-2 + e)(3c?3 - 1) + 3es?3(2c?2 " 1)) 7(2^3) eg , 
/iB^ = eci3 sin 26*12 61 — esi3 sin 26*12 eg . (68) 

The component B^ appears primarily due to Ami3, and the other component B^ vanishes if e = 0. 
Note that for two-flavors, or equivalently in the e = limit, B is completely contained in K*^^. 
Now, note the following structure in the equations of motion of a polarization vector: 

P^(u;, r) = H^(w, r) X P^(w, r) + H^(w, r) x P^(a;, r) , (69) 
P^(w, r) = H^(a;, r) x P^(a;, r) + H^(w, r) x P^(a;, r) . (70) 

It is clear from the above set of equations that if e = and one begins with P contained in K*^^, 
then P always remains in K^^, i.e. P^(a;,r) = identically. To investigate this case more closely, 
we write Eq. (j69|) for each component of P^ as ^ 

P3 = H4P5 - H5P4 , (71) 

P4 = H5P3 - H3P5 + \/3(H5P8 - HgPs) , (72) 

P5 = H3P4 - H4P3 + \^(H8P4 - H4P8) , (73) 

Pg = %/3(H4P5 - H5P4) . (74) 



^ In the following sections, the dependence of the Bloch vectors and the parameters on cj and r is implicit. 
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(75) 



Note that Pg = P3. This suggests that we could rotate our coordinates in the ea^eg plane by 
— 27r/3, so that Pg in the rotated frame becomes a constant of motion. While going to the rotated 
frame, the components X3 and Xg of any Bloch vector X transform as 

-1/2 -V3/2 
V3/2 -1/2 
The other components remain unchanged. 

This leads to the following simplified equations of motion for the two-flavor case: 

(76) 
(77) 
(78) 
(79) 

This is the two-flavor limit, where the state Vx does not participate in the evolution. This is a 
consequence of all the polarization vectors initially being contained in K^^. The rotated "tilde" 
frame can therefore be called as the "e — y" frame. 

The Eqs. ([76]) . ([77|) and ([78]) can be simply written as 

pes, ^ ^ey ^ pey ^ ^gQ) 

where the "x" can now be taken to be the usual cross product in a three-dimensional space 
spanned by {eg^, 64, e^}. This clearly exhibits the "gyration" of P about H, while the component 
of P along eg^ remains constant. The projection of P changes only along e*^ , which corresponds 
to t'e -H- flavor conversions. The problem is thus reduced to the two-flavor limit, for which 



P3 


= -2(H4P5 


- H5P4) 


P4 


= -2(H5P3 


- H3P5) 


P5 


= -2(H3P4 


- H4P3) 


Pg 


= . 
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analytical solutions have been discussed in literature 

In the two-flavor limit, it is observed that there are three qualitatively different kinds of motion of 
the polarization vector in the flavor space. The most familiar case is oscillations in vacuum/matter, 
where the neutrino-antineutrino density is small (/U <C oj) and each P((^) precesses about B with 
frequency cj. The other extreme is when the neutrino-antineutrino density is very large {fj, ^ uj). 
In such a situation, all P((^) remain tightly bound together and precess with the average uj of 
the ensemble, giving rise to synchronized oscillations. The intermediate regime (/u > lo) is when 
the P(w) remain bound together to a large extent, but have a tendency to relax to the state that 
has the lowest energy. The system is analogous to a pendulum/gyroscope that tries to relax to 
its vertically downward state, whatever state one might start in. This motion is called bipolar 
oscillation. 
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The motion changes quahtatively and quantitatively with the inclusion of the third flavor. There 
are two kinds of contribution due to the inclusion of the third flavor. First, we have some extra 
contributions to that depend on e, which changes the effective values of oj and ^13. These do 
not change the motion qualitatively. The second type of contribution is more interesting. It is due 
to the excursions of the polarization vectors into the Q subspace under the influence of B^. In 
particular, the length of is not preserved anymore. To see this clearly, we take the dot product 
of with Eq. ([69]) and that of P-^ with Eq. ^ to get 

|P^|2/2 = -|p^|2/2 = p^ .H^ X P^ . (81) 

We can clearly see that |P^|, which was a conserved quantity in the two-flavor case, no longer 
remains so. The non-conservation is proportional to |H^| and |P^|, both of which go to zero in 
the two-flavor limit. The addition of the third flavor makes the motion of the projection of P in 
the 63-63 plane fairly complicated in general, and we shall study it in some interesting regimes in 
Secimi 

F. The three-flavor solution 

In this section we extend the method presented in the last section to include the leading correc- 
tions due to the mixing of the third flavor. Let us illustrate our prescription in the vacuum limit, 
where the matter effects as well as the collective effects are neglected. The prescription will later 
be easily generalized to flnite matter densities and significant neutrino-neutrino interactions. 

Prom Eq. (152p . the Bloch vector B may be decomposed as 

wB = /iwB^^) + /iewB(2) + /iea;si3B(3) (82) 

with 

b(i) = sl^ 63 - 2(34 - 1)/(2V3) eg + (1 ~ esl^) sin 2^13 64 , (83) 

B^') = -(c?2 - c?3S?2) e3 + - 1)/(2V3) eg + C13 sin 2612 ei , (84) 

b(3) = 3si3(2c?2 - l)/(2\/3) eg - sin 2^12 eg . (85) 

Note that B(^) lies completely in K^y, B^^) in K'^^, and B(3) in K^?'. 

In Fig. [21 we show three coordinate frames e — x,e — y and x — y in the e3-eg plane. These 
frames are defined such that, if P is the projection of P in the es-eg plane, the components 
B(i)^ B(2)^ 3(3) ([85]) separately cause P to move along e^^,e^^,e^^ respectively. In order 
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FIG. 3: Useful coordinate choices on 63-63 plane 

to reduce the motions due to B*-"*^-*, B*^^-*, B^'^^ separately to two flavor problems as in Sec. Ill E] we 
write 

B^y = RB(i) , B^^ = B(2) , B'^y = R2b(3) , (86) 

where R is the rotation matrix in Eq. (j75p that rotates the X3 and Xg components of a Bloch 
vector in the 63-63 plane by — 27r/3. The vectors B'^^,B^^,B^^ are then simply B^-'^^ B(^\ B^^^ in 
the frames e — y,e — x,x — y respectively. We can then write Eq. (f82]) as 

(87) 



wB = uj^yR-^ B^y + w^^B'^^ + u'^yR-^ B'^y , 



with the "frequencies" defined as 



UJ 



ey 



and the "magnetic fields" as 



heoj 



xy 



heujsi3 sin 2^12 



B^y = cos 2^13 63 + (1 - esj^) sin 2^13 - (1 - 3sls)/{V3) eg 
B'" = -(c?2 - c?3S?2) 63 + ci3 sin 2^12 el" + (34 - 1)/(2V3) eg , 

^xy ^ _^xy _ ^s^3 cos 2^12/(4 sin 2^12) 63 . 



(88) 

(89) 
(90) 
(91) 



The vectors 64,61,66 are the directions transverse to the 63-6g plane that are relevant in the 
three frames, and can be written as e^^, e^f, 6^^ respectively. The -B°^ are normalized such that 
1 + 0{e, The separate motion due to each B^^ is then a precession about 



I Da/3|2 



\B 



a/3 1 2 



-S3 63 + e_|_ with a frequency w"'', where the half-angle of the cone is given by tan^"^ 
\Bf/Bf\. 
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The net motion of the polarization vector may be interpreted as the combination of two-flavor 
precessions about e^,e^ and e^^ respectively. It can immediately be seen from Eq. (|55j) that 

\oj^y\ > Iw^^l > Iw^^l , (92) 

i.e. the precession frequencies are hierarchical. Therefore, the motion due to slower frequencies 
may be neglected over short time scales. More precisely, if we coarse-grain the equation of motion 
Eq. (j58p in r over scales corresponding to uj^^ , the effects of oj^^ and uj^^ are negligible. The slowest 
variation in the solution is due to cj^^, which modulates the faster motion due to w'^^, which in 
turn modulates the motion at still shorter scales due to cj*^^. 

Let us denote the evolution of P(r) under the action of B*^^, B*^^, B^^ by the operators 
S'^^(r), S'^^(r), S^^(r) respectively. As long as the condition in Eq. (1921) is valid, we can write 

P(r) = S^y{r) S'=^(r) S^^(r) P(0) , (93) 

where the evolution operators are of the form 



1 



r/(w'=^,0^^,/i,r) 
1 



(95) 



S-(.)^R-J^(^"''J^'^'^)°|R^ (96) 

Here 77(0;"'', 0"'^, r) are the evolution functions that can be calculated in a two-flavor approxima- 
tion using the results in previous literature. In general, the frequencies of these evolution functions 
are determined by uj'^^s and the amplitudes are determined by the effective mixing angle ^"'^s. 
Each evolution operator S"'' takes the state P to the respective a — (3 frame in which Pg'^ stays 
constant and Pg^ undergoes precession, and brings P back to the 63-63 frame after precession. 
Note the matrices S"'^ are not unitary. The order in which they are operated should be such that 
the slower oscillations effectively act like an amplitude modulation for the faster oscillations. 
It is easy to calculate Pu^v^ using Eq. (j50|) as 



where P(r) is given by Eq. (I93p . and P3^,Pg^ are components along and 6g^ respectively. If 
neglect effects of the slowest frequency uo^^ , the expressions for Y'^ir) and Pg^(r) may be written 
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as 



P^^ (r) = r?(a;^^r^M,r) (^-^r/(.;-, r^ ^, r)P3(0) - ^PsWj , (98) 

PfW = (^+^r?(a;-,r^/z,r)P3(0)-^P8(0)^ . (99) 

In the presence of ordinary matter and when the collective effects may be neglected, the same 
prescription stays valid, simply with the replacements 

ei — )• cos ("(r) ei + sin ('(?') 62 , 64 — cos ({r) + sin ("(r) 65 (100) 

with ^(r) = V{r')dr' . It may be seen from Eq. (j6ip that these replacements take B to B*"*, 
so that the effect of MSW is taken into account by going to the interaction frame. As observed 
in Sec. IIIDl fast oscillations with a frequency ~ V{r) will average out the sinusoidal terms, thus 
decreasing the contribution from the transverse components of B^*"*)"^. 

When the collective effects dominate, since the collective potential Hyu{r) in Eq. ()36p is inde- 
pendent of energy, neutrinos of all energies precess with a common frequency in all the two-flavor 
subspaces. The motion is therefore similar to the vacuum case discussed above, with the replace- 
ment uj — )• {uj) as given in Sec. IIIIBI 

We have thus completed our program of expressing three-flavor effects purely in terms of two- 
flavor effects. The r-dependent functions r/(r) are known analytically for oscillations in vacuum 
and for synchronized oscillations, where we can explicitly check our ansatz. 

In the case of bipolar oscillations, the situation is more complicated since these are not sinusoidal 
oscillations, rather P remains almost static for a period of time and swings through the lowest 
energy state in a rapid burst. As a result, the fast- or slow-ness of bipolar oscillations as compared 
to the other precessions is time dependent. We therefore can obtain a qualitative understanding 
of bipolar oscillations in the three neutrino framework, but only a heuristic form of the analytic 
solution. 



III. FLAVOR CONVERSION MECHANISMS IN THREE-FLAVOR FORMALISM 

In this section we illustrate the three-flavor effects in some simple examples, where we take con- 
stant matter density and box-spectra for neutrinos and antineutrinos. We explain the three-flavor 
features therein analytically using the "es-eg" triangle diagrams. The insights gained thereby will 
allow us to understand the more complicated flavor conversions in realistic supernova simulations 
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in Sec. IIVI For our numerical evaluations in this section, we fix lAm^t^l = 2.5 x lO'^ eV^ and 
9i2 = 0.6. We also choose a box-spectrum for the the neutrino flux i.e. constant over the energy 
range E = (1-51) MeV, and zero elsewhere. 

A. Vacuum and MSW oscillations 

We start with looking at neutrino oscillations in vacuum/matter, with no collective effects. 
Although this situation has been analyzed in literature in great detail, we illustrate it here in order 
to familiarize the reader with the analysis in terms of F'^ and the "es-eg" triangle. This 
triangle, shown in Fig. HI helps in understanding the three-neutrino features of flavor conversions. 
The projection of P on the es-eg plane represents the flavor content, the allowed region being 
an equilateral triangle. The three vertices of the triangle represent the three states fe, I'x and i^y 
(anticlockwise, from top right). States that lie on the edges connecting them are admixtures of 
only those two flavors. The interior of the triangle represents states that are admixtures of all 
three flavors. Quantitatively, for any point on the triangle, the fraction of the neutrinos in flavor 
a is proportional to its distance from the edge opposite to the Va vertex, as shown in Eq. ()50p . 

In Fig. m we show the quantities P3^,Pg^ and Pu^ue as functions of the radial coordinate r. For 
illustration, we start with a pure z/e flavor, which corresponds to (P'^,P'^) = {—l,l/\/3). The 
following observations may be made from the figure. 

• The oscillation frequencies depend on the neutrino energy. However in the triangle diagram, 
the locus of P for all energies is identical for oscillations in vacuum (therefore, the thin and 
thick lines overlap). Different energies travel along this orbit at different, but constant speeds 
proportional to 1/E. In matter, the mixing angle begins to depend on the energy and thus 
the orbits are different for different energies. 

• The flavour evolution has two main frequency components. The fast oscillations with fre- 
quency uj = Am^^jjj/(2p) and the slower ones with frequency eu = AmQ/{2p). 

• If a; and eco were commensurate, the orbits in the triangle would be closed curves. However, 
that is a fine-tuned situation. In general, if e is not rational, the orbits do not close, but 
drift parallel to themselves periodically. Indeed, the orbits are analogous to the well-known 
Lissajous figures. 

• Pg^ only has slow modes corresponding to the frequency eoj. These slow oscillations modulate 
the amplitude of the upper envelope of jPg^j because the maximum value that jPg^j can take 
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FIG. 4: Neutrino oscillations for E — 20.0 MeV and 29.6 MeV (thin and thick lines respectively). To 
emphasize the nature the oscillations, we choose e — 1/5.1 and ^13 — 0.2. Oscillations in vacuum and 
matter are shown by dotted (blue) and undotted (red) lines respectively. For matter, we choose normal 
hierarchy and A = 0.3 km~^. 

is reduced when Pg^ deviates from its maximum value of l/VS- The above can be clearly 
seen from the triangle diagram. 

• Pg^ oscillations involve both frequencies, w and eoj. The maximum deviation of Pg^ from 
unity is governed by the amplitude of modulation of its upper envelop (which depends on 
sin^20i3) and the amplitude of faster oscillations superimposed on it (which depends on 
sin2 2012). 

• In the two-flavor limit we ignore the mixing with v^, and as a result Pg^ remains constant. 
In the triangle, this corresponds to the motion being confined to a line parallel to the 
axis. Indeed, the effect of the third flavor is to extend the motion of P to the entire triangle, 
as opposed to only along a line. The deviation of P from this line quantifies the extent of 
three-flavor effects. 

• The amplitude of oscillations can be read off from the triangle as the extent of the orbit 
along the Ve-i'y edge (2sin^20i3) and along the fg-z^x edge (2sin^20i2). 

• In the presence of matter, mixing angles are suppressed or enhanced depending on the energy 
and matter density. For A ~ euj, the MSW resonance occurs, and the effective mixing angle 
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becomes almost maximal, as it happens for the low energy mode shown in the figure. At 
X ^ eoj, the state Vx decouples because of the suppression of the mixing angle in matter, 
making this an effectively two-flavor ^ T^y problem. The oscillations in Pg^ have vanishing 
amplitude and the motion in the triangle is restricted to the Ve~i^y edge. 

• At even larger matter densities, A ^ w, the amplitude of Ve ^ '^y oscillations, which is the 
amplitude of Pg^ oscillations, starts decreasing and the motion in the triangle becomes more 
and more confined to be near the vertex as in the case of the high energy mode shown in 
the figure. 

All the above features may be understood analytically through Eqs. (j97|) - (j99]) and the two-fiavor 
evolution functions 



The above expressions are approximate, since we ignore the slowest frequency modes (depending 
on uj^y) and assume complete factorization. We find however, that these expressions agree with 
the numerical solution reasonably well . 

In the case of finite but constant matter density, we use the angles 6°^^ and frequencies cj"^ in 
matter, both of which are energy dependent. Note that the amplitudes in this case are proportional 
to 2 sin^ 20"^ in matter and can be maximal (spanning a full edge of the triangle) when there is 
an MSW resonance. 

When the matter density encountered by the neutrino varies such that neutrinos pass through 
an MSW resonance, they undergo fiavor transitions with adiabaticities depending on their energy, 
the relevant mixing angle and the matter profile. In the limit of a small mixing angle, a completely 
adiabatic H resonance is represented by a refiection of the neutrino state about eg^ in the 63- 
eg triangle. A non-adiabatic H resonance corresponds to a state that tries to move towards 
this reflected point, but does not completely succeed. Passage through the L resonance similarly 
corresponds to a reflection about Cg^. 




(101) 



(102) 
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FIG. 5: Synchronized oscillations for neutrinos oi E — 20.0 and MeV 29.6 MeV, which overlap completely. 
We choose e = 1/5 to emphasize the nature of oscillations, 613 = 0.2 and fj, = 100 km^^. Oscillations in 
vacuum and matter are shown by dotted (blue) and undotted (red) lines respectively. For matter, we choose 
normal hierarchy and A = 0.5 km^^. Note that the orbits on the triangle are the same for different energies. 

B. Synchronized oscillations 



At extremely large neutrino densities, it is expected that neutrinos of all energies oscillate 
synchronously with a common frequency (w"^) about B"'', given in the two- flavor case by 58] 

(-"^) = — , (103) 

where D's are the moments defined in Eq. (j63p . The frequency {uj"^) crucially depends on the 
neutrino energy spectrum. The box-spectrum that we have chosen corresponds to (o;*^^) ~ 0.49 
km^^. In Fig. Owe show Pg^, Pg^ and Pv^u^ as functions of the radial coordinate r for synchronized 
neutrino oscillations. 

The following observations may be made from the figure: 



The observations in Sec. IIII Al remain true, except that neutrinos of all energies oscillate with 
a common frequency in vacuum in the two flavor limits of each of the a — /3 subspaces. The 
response of all neutrinos to the neutrino-neutrino potential is thus identical. 

Even in the presence of matter, the synchronized oscillation amplitude is independent of 
energy, unlike what happens for non-collective oscillations. 
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• The amplitude of the slower oscillations is almost maximal because, in the chosen example, 
A ~ e(w). 

• The orbits drift periodically, even if oj and eu are commensurate, because {u) and (eco) 
are not commensurate in general. This is due to corrections to Eq. ()103|) arising out of 
incomplete synchronization. 

The above observations are explained analytically along the same lines as the vacuum / MSW 
case. The two-flavor evolution functions are given by 

7?K^r^0,r) = l-2sin2(2Msin2(^^^) , (104) 

.ex oex n „A _ 1 o„;„2/o/) \ „;„2 f he{uj)r 



r?(w'^^r^O,r) = l-2sin"(20i2)sin" . (105) 

In the plots we see that fast oscillations have wavelength 2ir/u} ~ 12 km. This matches the value 
of (w^^) calculated from Eq. (fT03|) . 

In the presence of a finite matter density, the MSW potential A also takes an effective average 
value given by 



58] 



(A) = A ^ . (106) 
Naturally, the mixing angle is also the same for all energies, since 

((A)/(a;"0 - cos 20"/3)2 + sin^ 20"/3 ' > 

Thus not only the frequency, but also the amplitude of oscillations is universal in the synchro- 
nized limit. The MSW resonance is collective, occuring with the same adiabaticity for all neutri- 
nos/antineutrinos at the same A when the relevant condition is met, as was shown in the two- flavor 



case 
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59|. The factorization shown in Sec. Ill El allows the result to be extended to the three- flavor 



situation. 



C. Bipolar Oscillations 

When the hierarchy is inverted and there are comparable numbers of neutrinos and antineutrinos 
in the system, i.e. /i|D| ~ i^|B|, the influence of the u and fx terms in the equations of motion 
depends crucially on the relative orientation of D,B and the magnitude of D itself. This subtle 
interplay gives rise to bipolar oscillations. 
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Many of the notions about bipolar oscillations in the two-flavor formalism [33, l33| remain valid 
with three flavors, since they do not depend on the number of flavors, or equivalently, on the 
dimensionality of the Bloch vectors. The system is best understood in terms of the "pendulum 
vector" Q defined in the interaction frame as 84] 



Q = S - -B , (108) 

in terms of which the equations of motion are 

Q = X Q - -B , (109) 
B = uBxQ. (110) 

The antisymmetry of the generalized cross product in Eq. (|57p implies that even in the case of 
three flavors, |Qp and D • B are conserved for large /x. 

In the two-flavor case, the motion can be understood in terms of a spherical pendulum 
with the total energy given by (^B • Q -|- /i|Dp/2. For normal hierarchy, the pendulum is stable 
and executes only small oscillations. For inverted hierarchy, however, the system behaves like 
an inverted pendulum, which tries to relax to its stable position. The polarization vectors then 
remain almost static, but periodically dip to the configuration with the lowest potential energy 
B • Q. Thus for inverted hierarchy, one can have a large flavor swap during the dip. The duration 
between successive dips is given by r''*^ « y^cJ/IlQI with logarithmic corrections depending on 9 
and A. Since fj, > uj, individual P remain bound to each other, and therefore behave identically to 

Q 

Addition of a third flavor may change the behaviour significantly, as we show in Figs. [6] and 
[7] for two extreme values of A. We consider the case of inverted hierarchy, and a box-spectrum 
of energies E = (1~51) MeV with the number of antineutrinos as (1 — a) times the number of 
neutrinos, with a = 0.2. Given that the hierarchy in the solar sector is normal, we expect bipolar 
effect only in the e — y subspace, combined with usual neutrino oscillations in the e — x subspace. 
The following observations may be made from the figures: 

• The evolution of both Po^ and Po^ consists of a series of bipolar "kinks" as in the two-flavor 



3 Clil^ X g 

case 



32l |. modulated by an envelope with the frequency {co'^^). The evolutions for neutrinos 
and antineutrinos closely follow one another, which is expected from the conservation of 
B D. 

Signiflcant three-flavor effects are present for small A, since the whole triangle is seen to 
be fllled with oscillations, forming a "petal structure" (Fig. [6]). It may be interpreted as a 
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FIG. 6: Bipolar oscillations at small A for neutrinos (dotted, blue) and 20% fewer antineutrinos (undotted, 
pink) of different energies, which almost overlap. We choose inverted hierarchy, |e| = 1/30, 6*13 — 0.01, = 
10 km^^ and A — 0.001 km^^. Note that the plots are the same for different energies, because of strong 
collective behaviour. 

combination of slow Vf. oscillations and bipolar oscillations that tend to take the state 
towards Vy in periodic bursts. 

• The extent of motion towards Vy depends mainly on the asymmetry a, whereas that towards 

depends on sin^ 2^12. 

• For large A (Fig. [7]), the oscillations in the e — x sector are suppressed since the effective 
mixing angle Oyi in matter becomes small. The amplitude of the bipolar motion is however 
not affected substantially. 

Bipolar oscillations (even in the two-flavor limit) do not have a sinusoidal form, hence they 
are not associated with a fixed frequency. They may be looked upon as a combination of a low 
frequency (during the time that the Vy component is stationary, which we shall call the A phase) 
and a high frequency (the sudden dip towards Vy^ which we shall call the B phase). Therefore, 
our prescription in Sec. Ill F] has to be applied with care. Note that the order of evolution matrices 
in Eq. ()96p is supposed to be in the decreasing order of frequencies. Even if we neglect the slow 
evolution due to B^^, strictly speaking during the A phase, one should use the order S^^S^^ and 
during the B phase, the order should be S^^S*^^. However, we find numerically that the evolution 
gexgey (j^Qgg^y matches the three-flavor solution over the complete evolution. This therefore may 
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FIG. 7: Bipolar oscillations at large A for neutrinos (dotted, blue) and antineutrinos (undotted, pink) of 
different energies, which almost overlap. We choose inverted hierarchy, |e| = 1/30, ^13 = O.OI,/i — 10 km~^ 
and A — 0.3 km^^. Note that the plots are the same for different energies, because of strong collective 
behaviour. 



be taken to be the heuristic solution for the bipolar oscillations in the three-flavor case. 

We have not considered normal hierarchy, in which we expect that starting with we'll have a 
stable system that will not undergo bipolar oscillations, whereas starting with Vx or i/y, we'll have 
independent bipolar oscillations towards I'e- It will be interesting to analyze the details of such a 
scenario, however it is beyond the scope of this paper. 



D. Spectral splitting 

As a system of neutrinos and antineutrinos transits from the collective regime {fi ^ a;) to 
vacuum ~ 0), the polarization vectors P keep trying to align with H in the adiabatic approxi- 
mation. Due to the conservation of B • D, as shown in Sec. IIII C| this alignment is not possible for 
all P. Indeed, neutrinos with high energy need to flip over and anti-align with B (which equals H 



m vacuum) 



341, 



35|. This leads to a sharp split in the energy spectrum, with the high energy 
getting completely converted to the non-electron flavor and vice versa. 

A crucial requirement for the splits to develop is the preparation of the sytem for the split by 
the generation of components of P that are transverse to B. Bipolar oscillations do this easily for 
inverted hierarchy, independent of matter effects. For normal hierarchy, in the presence of large 
matter effects the oscillations are suppressed, but an MSW resonance can prepare the transverse 
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FIG. 8: Spectral splits at large A for neutrinos (dotted) and 33% fewer antineutrinos (undotted) with a box- 
spectrum of energies E = (1-51) MeV. In the Pi/^i/^ plot, the energy of neutrinos (antineutrinos) increases 
(decreases) top downwards. The energies (in MeV) of the modes, shown in the figure, are 1.0 (Red), 1.5 
(Blue), 3.5 (Green), 12.5 (Pink) and 32.0 (Grey). We take inverted hierarchy, |e| = 1/30, 6*13 = 0.01, 
fi = 10^ (50/r(km))^ km^^ and A = 10 km^^. In the e^-eg triangle, the evolution is always along the i^e-i^y 
edge. 

components. 

For illustration, we choose two situations, with large and small A (Fig. [8] and [9] respectively) 
and the hierarchy is taken to be inverted. We choose the box-spectrum for and Ue energies and 
the flux asymmetry a = 0.33. We observe the following from the figures: 

• For large A, there is only a single split for neutrinos, which can be seen in Pg^. The split is 
not visible in the triangle since the neutrinos are confined to the Ve^Vy edge. However, the 
low energy neutrinos move towards Vf, and the high energy ones towards Vy. 

• for small A, the split is not only in Pg^ but also in Pg^. There also are oscillations with large 
amplitudes. Some neutrino states drift towards and ultimately reach Vy^ while the others 
keep oscillating between and v^- 

The above observations can be understood as follows. For large A, the solar mixing angle is 
suppressed and hence the problem reduces essentially to a two-fiavor one in the e — y subspace. 
Thus, the split is only in Pg^. The split happens in neutrinos since there are more neutrinos than 
antineutrinos at any given energy. For small A, in addition to the above split, there are large 
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FIG. 9: Spectral splits at small A for neutrinos (dotted) and 33% fewer antineutrinos (undotted) with a 
box-spectrum of energies {E — 1-51) MeV. The conventions for lines is the same as that in Fig. [5] We 
take inverted hierarchy, e — 1/30, 613 ~ 0.01 /i = 10^ (50/r(km))'* km^^ and A = 0.1 km^^. In the es-eg 
triangle, we show only some of the representative energies that have different behaviours. 

Vf. o oscillations, which give rise to a split that is observable also in Pg^, which was absent for 
large A. 

A detailed understanding of the spectral splits in the three-flavor case, including predictions for 
the positions of the spectral splits has been obtained 

IV. COLLECTIVE EFFECTS ON NEUTRINOS IN SN 

Collective effects are likely to be important in the context of neutrinos emitted from a SN. The 
number density of neutrinos and antineutrinos streaming off the neutrinosphere is quite large, so 
dominates over A and cj upto a radius of a few ten or hundred kilometers respectively. Therefore, 
it is likely that one or more flavor conversion mechanisms identified in Sec. IIIII come into play in 
the different regions inside the star. 

In this section, we study the effect of collective oscillations and their interplay with subsequent 
MSW transitions inside a SN. To illustrate the nature of these effects, we numerically solve for the 
flavor evolution equations by taking a realistic SN density profile, including the collective effects. 
We present results primarily for inverted hierarchy, because collective effects are not expected to 
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play a significant part for normal hierarchy ^. In the numerical study, we take |Amf3| = 2.5 x 10~^ 
eV, |e| = 1/30, 9i2 = 0.6, and two representative values of 6*13, viz. 6*13 = 0.001 and 0.1. We 
then apply the formalism developed in Sec. |TI] to the case of neutrinos streaming from the SN 
neutrinosphere and identify the regimes where different flavor conversion mechanisms are at work. 
This allows us to explain the features in the observable neutrino and antineutrino spectra, and 
understand the three-flavor effects. 



The SN model is defined by the following choice for the emission geometry, initial flavor de- 
pendent spectra and fluxes, the collective potential and the matter density profile. We would like 
to emphasize that these choices are canonical and more importantly, the specific value of the lu- 
minosity or the spatial dependence of the collective potential does not affect results significantly. 
Any large initial value of // (such that it exceeds to) and its slow decrease with r gives almost 
identical results. In other words, the results are not fine-tuned. However, the results would depend 
on the flavor spectra and the matter density profile, as these determine the initial conditions, the 
collective potential and the effective mixing parameters. 



Neutrinos with different energies and flavors start freestreaming at different r, but flavor evo- 
lution does not start until much later. Thus the radius of the neutrinosphere rg is used only to set 
the initial conditions. We therefore use the "bulb-model" of neutrino emission from the SN as dis- 



cussed in ^] with a nominal neutrinosphere at tq = 10 km. We assume steady-state half-isotropic 
emission from the neutrinosphere. 



A. SN model and parameters for numerical simulation 



1. Emission geometry 



n 



2. Initial spectra and fluxes 



The flavor-dependent primary neutrino spectra at rg are parametrized as [65 1 




exp -{^a + 1) 



E 



(111) 



^ However, during the neutronization burst phase of an 0-Ne-Mg SN, this need not be true. 
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where N{^) = (l+^)^"'"^/r(H-^). This spectrum is normahzed such that dE^^ Pua{Eua) = 
and has the average energy {Ey^). The above parametrization has the advantage that the spectra 
can be analytically integrated, including the effects of spectral pinching through The number 
flux is given by = Ly^j {Ey^ ^ where Ly^ is the luminosity in the flavor v^. We remark 
that an equipartition in the luminosity is often assumed for simplicity. For illustration, we choose 
the above parameters as 

Ly^ = 1.5 X lO^^ergs/sec , = 3 
[Ey^^ = 10 MeV , (^p,) = 15 MeV , (^,,,^, p,,^) = 20 MeV . (112) 

Remembering that ~ p = |Amf3|/(2a;), we can rewrite the above information in terms of w, 
if desired. Combining Eq. (I48p with the definitions of moments in (I63p . allows us to calculate the 
values of D(ro), S(ro) and D(^)(ro) for the above spectrum as 

- {EyME...) + {E,^) (MEyJ + (EyJ) " 1 1 ^^^^^ 

^( . _ {{EyJ + {EyJ){EyJ-2{EyJ{EyJ ^ _ 2 

^ " + {E,J) {EyJ + MEyJ{E,J - 11 ^''^^ 

X 2Aml. l/{Eyf + l/(Ey)^ -2/(Ey )^ -, , , 

ro = ^ '\, '' , , /,L ee = 0.215 ee km~^ 115 

^ ' 3 l/{EyJ + l/{Ey^)+A/{EyJ ^ ^ 

Using the above expressions, (w) = D • D(^)/|D|2 is calculated to be 

(w) = 2.37 km-i , (116) 
which allows us to write {lo'^^) = {oj) and (lo^^) = e{uj) in terms of (w), as per Eq. ()103p . 

3. Collective potential and matter density profile 

The collective potential for r > tq for the choice of parameters in Eq. (I112p is given by 

Ai(r) = 0.45 X 10^ g{r) km"^ , (117) 

where g{r) is given in Eq. (I28p . For illustration, we choose the shock- wave simulation inspired 
density profile that corresponds to ^ 

A(r) = 1.84 X lOV?^^'"^ km~^ . (118) 

The profiles of A(r) and //(r) are shown in Fig. [TOl 



This is the same as the one used in [l^ at t=4 sec. 
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FIG. 10: The profiles of A(r) and ^{r) for tlie SN model chosen in this section, along with the bands for 
the MSW resonances H and L. Also indicated are the terminal values of r where synchronized /bipolar 
oscillations for the e — y and e — x flavors take place. 

B. Flavor conversions inside a supernova 



In this section we solve for the evolution of the neutrino density matrix numerically for the chosen 
density profile, and show the neutrino flavor conversions. We expect synchronized oscillations in 
the region where ^ > A{u)'^y) Sg^ / (Dg^)^ 



208 km"^ 

in our example. In the further region till ^ ~ {uj^^) / Dg^ 



48 1, which corresponds to rgyn ~ 30 km 



26 km" 



481], which corresponds to 



'"'hip ^ i^y bipolar oscillations are expected. Beyond this region the spectral split in 

the e — y sector should develop, and subsequently MSW resonances should take place. Similarly 



we calculate for the e — x flavors, the relevant values of r 



syn 



68 km and r^^^ 



114 km for 



approximate boundaries of synchronized and bipolar oscillations in the e — x sector. However, no 
bipolar oscillations take place in the e — x sector since the corresponding hierarchy is normal. In 



Fig. [Tot we show the positions corresponding to r 



syn ! ''"trip ! ''"tyn 



and rf? , 



Small 9i3 



Fig. 111! shows the flavor evolutions in terms of P3^,Pg^ and the 63-65 triangle for neutrinos as 
well as antineutrinos, for ^13 = 0.001. This small value of ^13 ensures that the MSW resonance 
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FIG. 11: The flavor evolution of representative energy modes of I'e (dotted) and (undotted) for the density 
profile in Fig.[TUl with 6'i3 = 0.001. In the p^^v^ plot, the energy of neutrinos as well as antineutrinos increases 
top downwards. The energies (in MeV) of the modes, shown in the figure, are 2.5 (Red), 3.6 (Blue), 9.4 
(Green), 13.3 (Pink) and 50.0 (Grey). In the triangle plot, the bold line passing through is where all the 
neutrino and antineutrino states initially lie. 

H in antineutrinos is nonadiabatic, so that the effects of this resonance are not felt. One can 
then cleanly identify the collective effects. We make the following observations and interpretations 
based on the figure: 

• All the neutrinos and antineutrinos initially lie on a line passing through Uf. in the es-eg 
triangle. This is because the initial conditions are taken to be symmetric in and Uy. 

• The flavor evolution starts only at r = 40 km, which is slightly beyond rfy„. Before this 
point, the oscillations are synchronized, with a vanishing amplitude since A ^ {uj^y). 

• Between r = 40 and 60 km, V(, o Vy bipolar oscillations are observed as rapid dips in Pg^, 
and consequently in pv^y^. These oscillations vanish when r > r^?^. 

• Around r ~ 60 km, a spectral split develops in neutrinos along Pg^. The spectral split tends 
to keep the low energy neutrinos at their original position, while taking the high energy 
neutrinos as well as almost all antineutrinos towards Pg^ = 0. ^ 



There seems to be a spectral split in antineutrinos as well, at very low energies. This is similar to the observation 
in [4^ . and may be the effect of nonadiabaticity in the splitting. 
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FIG. 12: Neutrino and antineutrino spectra at r = 5000 km for ^13 = 0.001. The e, x and y flavors are 
shown in red(soUd), green(dashes) and blue(dots). The thin hnes/dashes/dots are for initial spectra and 
thick ones for the final spectra. The and Vy spectra get swapped for > 7 MeV, whereas the lower 
energy j/g spectrum partially mixes with Vx- In the antineutrino sector, the and Vy spectra are almost 
completely swapped, while the Vx spectrum remains unaffected. 



• Between r ~ 100 — 1000 km, antineutrinos of different energies undergo the H resonance. 
However the resonance is highly nonadiabatic and does not cause any flavor conversion. 

• At r = 1000 km and beyond, the effects of the MSW resonance L come into play, resulting 
in Vf> -f-T- Vx conversion. Since the high energy neutrinos are already close to = 0, there 
is effectively no flavor conversion. However the low energy neutrinos tend to convert to Vx, 
which is seen as a movement parallel to the i'e~i^x edge in the 63-63 triangle. 

• Since all the flavor conversions can be understood as a net effect of two-flavor phenomena 
taking place in well-separated regions in the star, the flavor transitions in the 63-63 triangle 
are always parallel to the i'e~i^x edge or Ver^y edge. 

Thus, for a small ^13, the collective effects can be clearly identified, whereas the effects due to 
the H resonance are absent. We calculate the flavor evolution till r = 5000 km. The collective 
effects have almost vanished by this time. Further MSW resonances due to the shock wave 
as well as possible effects of stochastic density fluctuations or turbulence 



1J2O 



11 



2C] will govern flavor 

conversions here onwards. Our calculations thus provide initial conditions for neutrino spectra at 
this point. 

In Fig. [T2I we show the neutrino and antineutrino spectra at r = 5000 km. We see that 
with E MeV convert almost completely to Vy due to the spectral split, whereas lower energy 
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FIG. 13: The flavor evolution of same representative energy modes of j/g (dotted) and (undotted) for the 
density profile in Fig. [TUl with (?i3 = 0.1. The convention for the lines is the same as in Fig. [11] 

Vf. convert partially to at the L resonance. In the antineutrino sector, the and Vy spectra are 
almost all completely swapped due to the spectral split, while the Vx spectrum remains unaffected. 

2. Large 6*13 

At large ^13 values, the H resonance at r ~ 100-1000 km is adiabatic, and causes significant 
flavor conversions in antineutrinos. In Fig. [131 we show the flavor evolution for ^13 = 0.1. While the 
signatures of synchronized and bipolar oscillations as well as the spectral split remain identical to 
the ^13 = 0.001 case, the K resonance can be seen to change the antineutrino picture substantially. 
The conversions in the neutrino sector, on the other hand, are identical to the small Q\z case. The 
following observations can be made from the figure. 

• The spectral split gives rise to a complete Ve~'^y conversion, which takes antineutrinos to 
Pf = 0. 

• The H resonance again swaps the Ve^Vy spectra, thus undoing the earlier effect of the spectral 
split. This takes the antineutrinos back to their starting position in the triangle. 

• Antineutrinos are now not on the = line as in the small ^13 case. As a result, the 
large value of Q\2 causes substantial Ve^^x conversion as the neutrinos emerge from the L 
resonance region. 
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FIG. 14: Neutrino and antineutrino spectra at r = 5000 km for 613 — 0.1. The e, x and r/ flavors are shown 
in red(soHd), green(dashes) and blue(dots). The thin hnes/dashes/dots are for initial spectra and thick 
ones for the final spectra. The i^e and i^y spectra get swapped for E > 7 MeV, whereas the lower energy 
spectrum partially mixes with v^- In the antineutrino sector, the De and spectra are partially mixed, 
while the Dy spectrum remains unaffected. 



The neutrino and antineutrino spectra at r = 5000 km are shown in Fig. [T4l We see that the 
neutrino spectra have the same characteristics as for small ^13. In the antineutrino sector, complete 
i>e-Vy spectral split and the reconversion at the H resonance cancel each other, whereas the large 
value of 9i2 partially mixes the i^e-^'x spectra. 

The value of ^13 thus affects the De spectra substantially. At larger ^13 values, where the H 
resonance is more adiabatic, the 9^ spectrum is softer. The spectrum is also affected at large 
^13, as opposed to the small ^13 case. 



3. Summarized results 



It is thus clear that the neutrino fluxes that reach Earth from a SN, are very different from 
the primary fluxes. In particular for inverted hierarchy, we learn that the v^. and Vy spectra are 
exchanged above a certain split-energy Ec due to collective effects. For antineutrinos the swap 
occurs over the complete spectrum. In the normal hierarchy, collective effects do not have a 
significant effect. The MSW conversions cause further flavor conversions, and while the conversion 
probabilities have not changed from the traditional expectation, the primary fluxes entering the 
resonances are now vastly different. This leads to different flavor composition of the fluxes of 
neutrinos and antineutrinos arriving on Earth, than was traditionally expected. These fluxes can 
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be calculated using our understanding of collective effects and the level-crossing diagrams. At the 
detectors on Earth one is typically sensitive to the fg and/or Vf, flux, and so we summarize the 
expectations for these fluxes in Table HI The expressions in the table are able to describe all the 
features of fg and v^. spectra in Figs. [T2] and [HI 



Normal hierarchy 




Inverted hierarchy 




Ft = s?2 {PhF,^ + (1 - Ph)F,,) 


+ cl^F,^ 


po,s^\'l2F.^+'^l2F.. 


{E < E,) 
[E > E,) 


Ft = CI2F,. + s^Fv^ 




Ft^42F,^+c',2{(^- 


PH)Fi,^ + PhF^^) 



TABLE I: Neutrino and antineutrino fluxes arriving on Earth from a SN. 



We have taken the L resonance to be adiabatic. In the case of multiple H resonances, as may 
occur during the shock wave propagation or turbulence, Ph may be taken to be the effective jump 
probability (it may have a nontrivial dependence on energy and time). Note that Earth matter 
effects are present only when Ftp^ is a nontrivial combination of Fy^jp^ and Fy^^y^jy^^p^. 



V. CONCLUSIONS 



We have developed a formalism to analyze neutrino flavor conversion effects in the full three- 
flavor framework. It employs the Bloch vector representation for 3x3 density matrices, and 
naturally generalizes the spin-precession analogy to three flavors. In particular, it is capable of 
describing three-flavor collective neutrino conversion effects inside a core collapse supernova, like 
synchronized oscillations, bipolar oscillations and spectral split, which have till now been analyti- 
cally studied mostly in the two- flavor limit. 

We explicitly extend the earlier two-flavor analysis of neutrino flavor conversions inside the SN, 
which includes neutrino-neutrino interactions, to three flavors, where we neglect the CP violation 
in the neutrino sector. We use the modified flavor basis (z/g, i^^, fj^), which is rotated from the 
flavor basis (i/g, t'^, i^t) so as to get rid of the mixing angle ^23- We also work in the steady state 
approximation so that there is no explicit time dependence in the density matrix, assume spherical 
symmetry and half-isotropic neutrino source, and employ the single-angle approximation that has 
been shown to be valid in the two-flavor case. This leads to the equations of a gyroscope in eight 
dimensions, similar to the three dimensional gyrosope equations in the two-flavor case. 

In the three-flavor formalism, the density matrix is represented by an eight-dimensional Bloch 
vector P. However, the flavor content is determined only by the two components P3 and Pg 
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of P after evolution. Motivated by this observation, we propose the "ea-es" triangle diagram 
to represent the flavor content of any neutrino state by the projection of P on the 63-63 plane, 
which we have termed P. This not only allows us to visualize the three-flavor transformations in a 
convenient way, but also allows us to quantify the extent of three-flavor effects over and above the 
two-flavor results. 

A "heavy-light" factorization holds in the three-flavor treatment for certain initial conditions, so 
that the three-flavor results may be understood as the two-flavor results with Am^^j^^ modified with 
terms that depend on Amg. Indeed, in certain situations, the three-flavor neutrino conversions 
may be f actor ized into three two-flavor oscillations with hierarchical frequencies. In such cases, the 
three-flavor conversion probabilities may be constructed from two-flavor results by considering the 
modulation of higher frequency modes by lower frequency modes. 

We have compared our analytic results with the numerical ones for simple cases of an initial pure 
z/g state, constant matter densities and no collective effects, as well as for synchronized oscillations, 
and have found a good agreement even when we ignore the modulation due to the lowest frequency. 
The additional effect of the third neutrino in these cases is limited to the excursions of the orbit of 
P towards v^-, and has already been well studied (though without the Bloch vector treatment). In 
the absence of collective effects, though the evolution of all energies is different, the orbit of P can 
be seen to be an energy-independent quantity. In the synchronized case, neutrinos of all energies 
are seen to oscillate with a common frequency, and even undergo MSW resonances at the same 
matter density and with the same adiabaticity. 

In the case of bipolar oscillations, the addition of the third neutrino changes the situation 
significantly. The analytical results are not so easy to obtain, however the numerical results for 
an inverted hierarchy show a "petal" pattern in the 63-63 triangle, which can be explained by the 
combination of i/g ^ bipolar oscillations and i^x sinusoidal oscillations. The value of the 

MSW potential also plays an important role in determining the extent of the effect of the third 
flavor. This needs to be explored in more detail. 

The spectral split occurs in neutrinos in the inverted hierarchy when one starts with f^, owing 
to the unstable position of the eight-dimensional gyroscope in this case. The z^g above a certain 
energy, and almost all v'^, completely convert to I'y and I'y respectively. There arc no additional 
split effects from the introduction of the third flavor since the hierarchy in the solar sector is normal. 
This, however, could change if neutrinos are not in a pure state as they enter the bipolar region. 

We have simulated the neutrino flavor conversions numerically by taking a realistic density 
profile for the SN, and have shown the flavor conversions for inverted hierarchy and two ^13 values 
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in the cooling phase. In such a scenario, it is easily possible to identify regions where different 
collective as well as MSW effects dominate. We are able to predict the regions in which these 
effects take place, and our three-flavor formalism can explain the features of flavor conversions 
therein. We also point out an interplay between the collective and MSW effects. For example, 
the H resonance cancels the effect of the spectral split for antineutrinos, whereas the spectral split 
makes the L resonance irrelevant for neutrinos above the split energy. If the hierarchy were normal, 
the collective effects would be effectively absent in the cooling phase. 

In conclusion, a complete understanding of the neutrino flavor conversions inside a SN requires 
a three-flavor treatment. In this paper, we have developed a formalism to handle this analytically, 
provided a method to estimate the three-flavor probabilities using the two-flavor results in certain 
situations, and have pointed out an interplay between collective and MSW effects, which can be 
easily understood with the formalism. 
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Appendix A: Notation 

This paper uses 3-vectors, 3x3 matrices, 8-vectors, their components and their projections in 
two-dimensional planes. We have tried to be consistent in the use of fonts for all these objects. 
The convention followed for indices is: 

• k G {1, 2, 3} index mass eigenstates 

• a, P E {e, X, y} index modified flavor eigenstates 

• a,b,c & {1, 2, 3, 4, 5, 6, 7, 8} are SU(3) indices. 
The convention followed for symbols is: 

• Three-dimensional vectors are denoted by smallcase boldfaced letters, e.g. p,q, r,v. 
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• 3x3 matrices are denoted in the regular math mode, e.g. p, H, I, M,0, Rij ,U,V, Aa- Note 
however the exceptions, energy E and the Fermi constant Gp, that are also respresented in 
this font. 

• Sets of 3 X 3 matrices are denoted by blackboard-bold font, e.g. K, Q. 

• Eight-dimensional vectors are written in bold capital letters, e.g. B,D,P,S,X. Exceptions 
are Gq, which are vectors along the coordinate axes, and therefore conventionally written in 
smallcase. 

• The components of an eight-vector X are written as Xq. 

• A is a eight-dimensional vector whose components are 3x3 matrices A^. 

• Projections of eight-vectors on the ea-eg plane, e.g. P, as well as the rotation and evolution 
matrices on that plane, e.g. R, S are written in sans-serif font. 
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